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Using conformal field theory, we perform a complete analysis of the chiral six-point correlation 
function 

C{z) = {(f>l,24>l,2^1/2,o{z, z)4>l,2(j}l,2) , 

with the four 0i,2 operators at the corners of an arbitrary rectangle, and the point z = x + iy in the 
interior. We calculate this for arbitrary central charge (equivalently, SLE parameter k > 0). C is of 
physical interest because for percolation (k — 6) and many other two-dimensional critical points, it 
specifies the density at z of critical clusters conditioned to touch either or both vertical ends of the 
rectangle, with these ends 'wired', i.e. constrained to be in a single cluster, and the horizontal ends 
free. 

The correlation function may be written as the product of an algebraic prefactor / and a conformal 
block G, where / = f{x, y, m), with m a cross-ratio specified by the corners (m determines the aspect 
ratio of the rectangle). By appropriate choice of / and using coordinates that respect the symmetry 
of the problem, the conformal block G is found to be independent of either y or x, and given by an 
Appell function. 
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I. INTRODUCTION 

The methods of conformal field theory (CFT) [1, 2] allow calculation of the correlation functions of a variety of 
operators, which may, in many cases, be interpreted as physical quantities in the continuum limit of two dimensional 
critical systems. As the number of operators grows, the calculation becomes progressively more difficult, however, 
and in practice there are very few results for correlators with more than four operators. 

In this paper we present a full calculation of the six-point correlator 

Ciz) = {<PUmi2mi/2,0{z, -z)c^l2{RWl,2{R + i))7J , (1) 

in the rectangular geometry TZ := {z = x + \y £ C\Q < x < R,Q < y < 1}. The conformal dimensions used through 
out this article and the associated central charge are 

6-^ - (8-«)(3^-8) 
ni,2 = — — 'ii/2.0 - "-1/2,0 — 71 l^J 

8 -At (3k-8)(6-k) 
"1,3 = c= , 

where k is the Schramm-Loewner Evolution (SLE) parameter. Because of their positions the (pl 2 corner operators 
have an effective dimension hi 2 = 2/ii,2, i-c. twice the usual value (see subsection II B for more details on this point). 

C(z) is of interest because it may be used in a variety of physical models to determine the density of clusters 
attached to one, or both, of two distinct boundary intervals, when these intervals are each 'wired', i.e. constrained 
to belong to a single cluster. A recent paper [3] presents a calculation of C for percolation (k = 6) in a semi-infinite 
rectangle, and also in an arbitrary rectangle by assuming a certain y-independence specific to the rectangle. 

Here, we calculate C completely and without assumptions, for arbitrary k > 0. This is done by solving the 
differential equations that C satisfies. The main new step is a certain choice of coordinates. This choice allows us 
to derive the curious y-independence mentioned, to reduce the number of variables from three to two, and write the 
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solutions explicitly in terms of Appell functions. The application of these results to a variety of physical systems is 
considered in [4]. 

To begin, we express C in the form 

C{z)^f{z)F{z), (3) 

where f{z) is an appropriately chosen algebraic prefactor. With this choice it transpires that for a given rectangle, all 
possible solutions for F{z) can be written in the form F{z) = G{x) + G(y), where G is determined by G. Then, using 
certain elliptic functions of x and y as intermediate variables, we find an algebraic expression for /, and algebraic 
factors times Appell hypergeometric functions for G. The full expressions for / and G depend on the aspect ratio of 
the rectangle as well as x and/or y. 

The prefactor / is independent of the details of the physical system, while F changes depending on the particular 
observable associated with C{z). For a conformal block of (1), F{z) only depends on x or y. This surprising feature, 
originally observed numerically (see [3]), indicates the presence of some unknown symmetry. It also implies that in a 
given rectangle appropriate ratios of two C{z) with different physical meanings are completely independent of either 
X or J/, since the prefactor / cancels out. 

In a companion paper [4] we apply these results in several ways. First, to find the cluster densities for a range of 
critical 0(n) loop models, in both dense and dilute phases and equivalently, for critical Q-state Potts models, probing 
either FK or spin clusters. Second, we extend previous results for the factorization of correlations for percolation, 
described in [3] , to the critical models mentioned. Finally, for percolation, the density of horizontal crossing clusters 
in a rectangle with open boundary conditions on all edges is determined. 

In this work, section II calculates the correlation function (1) by solving the associated PDEs, which with proper 
choice of co-ordinates and prefactor / reduce to the Appell equations. Subsection II A gives the PDEs, choses / and 
coordinates, and presents the solutions for G, all of which are single conformal blocks. The interesting independence 
of G from one coordinate mentioned (see (15)) appears here. Subsection II B computes the form of the correlation 
function prefactor / in the rectangle, a not completely trivial task. Section III contains a summary of our results and 
some discussion. 

In Appendix A, we derive a relation between Appell functions that is useful in simplifying one of our formulas. 
Appendix B examines conditions for a set of conformal blocks to have a common y-dcpcndence in the rectangle. 



II. THEORY 



A. Solving the differential equations 



In this subsection we determine and solve the differential equations for the correlation function (1). The main new 
step in solving the differential equations is a certain choice of coordinates, given below. This choice allows us to derive 
the interesting y-independence mentioned, and to write the solutions in terms of Appell functions. 

To begin, we consider C in the upper half plane H := {w = u + it; | u e K, d > 0} 



C{w) = (01,2("l)(/'l,2(w2)$l/2,o(w, M')01,2("3)</'l,2(w4))l 



(4) 



using the methods of boundary CFT, then transform into the rectangle TZ. In H we can decompose $i/2,o(''^j ^) voio 
chiral components ^i/2.Q{'w)^i/2fi{w)- Then, by conformal symmetry one may write 



G{w) 



\w ■ 



/ll,3 — 2/li/2,0 



\{w - Ui){w - U2){w - U3){w - U4) 



{W - Mi)m43 (W - Mi)m43 ^21^43 
3'*i.3 \M3l(u4 — w) ' U3l(u4 — w) ' M31M42 



-F 



(5) 



where 



— Uj. In the next section we give our reason for choosing this particular form for C{w). 



Using standard CFT methods [1] F may be determined via the differential equations arising from the null state 
[3L_i^ — 2(1 + 2/11.2)^-2] 101.2) associated with each 0i^2 in (5). The presence of 4>i.2{ui), for instance, means that 
(5) is annihilated by the operator 

2hi/2,oRc[{w -uif] 2Re[(w-ui)9„] /ii,2 



\w — Ui\ 



\w 



Ul 



+ 



i=2 



Ml) 



2(1 + 2/11,2) 



(6) 



Next we let {wi, U2, U3, 1x4} — > {0,m, l,oo}, which means that under the conformal map to the rectangle (8) m is 
the cross-ratio m = (U12U34/U13U24) of the image points of the corners of the rectangle. (Note that m differs from the 
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standard modular lambda parameter, which is 1 — to here.) Thus one arrives at a differential equation for F{w, w, to), 
8(6 - k)ww + (8 - k)to(4(2to - (w + w)) + (8 - k)to(1 - w){l - w)) 







l2 

(1 — to)((8 — K)m{w + w) + 2'ww{Km — 4)) 
2Kmww 



dmF 



{I — w){{8 — k)w ~ Kw + 2kww) {1 — w){{8 ~ k)w — kw + 2kww) , . 

H n : H ; OtB-f* (7) 

Zkww 2kww 

- (1 - mfdmdraF - 2(1 - w){l - m)9^9,„F - 2(1 - m)(l - w)d^dmF 

- (1 - u')2a„a^i? - 2(1 - w)(l - w)9a9„F - (1 - wfd^d^F . 

We next transform (7) into rectangular coordinates via the conformal mapping 

^(z, to) = TO sn (z i^T'lm)^ , w{z^m) = m?,n{z K'\rrif' , (8) 

where K' := K{1 — to), with K := K{m) the complete elliptic integral of the first kind, and sn(-|TO) the Jacobi elliptic 
function with elliptic parameter to. The factor K' appears because our TZ has fixed height of 1, which differs from 
the standard rectangle used to define the Jacobi elliptic functions. The aspect ratio is given by 

K(m) , , 

a'(to) 

which is the inverse of the standard elliptic aspect ratio. Conversely, the aspect ratio R specifies the elliptic parameter 
TO via 



(0,e- 



(0,e- 



TrR\ 

(10) 



The transformation (8) maps the corners of the rectangle TZ, starting at the origin and proceeding counterclockwise, 
into the w-plane points (0, 0), (to, 0), (1, 0) and (oo), respectively. 
We now introduce the real coordinates 

= sn{xK'\m)'^ , and ^ ~ sn {y K' \1 ~ m)'^ . (11) 

This choice of co-ordinates is a key step, as we will see. It simplifies the equations, leading to a solution of the PDEs. 
In addition, our results are either algebraic or hypergeometric when written in terms of ^ and ip. 
The conformal transformation now becomes 

^(1 - (1 - mjjj) - (1 - 0(1 - mQijl - ^) , . 2Ve(l - 0(1 - ^0^(1 - - (1 - m)ij) 
w[C,V,m) ^ m — -— hiTO — — ^ . 12 

The coordinate to within a factor to, is the half-plane image of the projection z ^ x onto the bottom edge of 
the rectangle, i.e. w(x) = so that ^ S (0, 1). The coordinate "0 is slightly more complicated and is determined by 
first taking the projection z i— ^ iy on the left end of the rectangle. The half-plane image of this point is in the interval 
(— oo, 0), and we define ip via w{iy) =: m-tp/lil) — 1) so that -0 G (0, 1). We will see in the following paragraphs that 
this choice of coordinates introduces a useful ^ -fT- -0 symmetry and allows us to compare directly with results from 
[3]. 

Transforming (7) into these coordinates (with the help of Mathematica) yields 

= (8(6 - K)(e + -0 - ^0)^ + (8 - Kf{l - TO$ - + mi)f (13) 
+8(8 - n) ((1 - 0(1 - ^) - ^0(1 - mO(l - ^' + W^))) F 
+ 2kto(1 - to) (((8 - K)(e - V') - 8^0(1 - 2to))(2 - ^ - + ^1/^) - (1 - 2to)k (^^ - + V-^ - C0^)) S^F 
+ 4«;01 - (4 (1 - + to2^0) - (k - 2< + 40(1 - V' + mi]:f - (k - 4)^to(1 - -0 + mtp)) d^F 
+ 4k0(1 - 0) (4 (1 - C + (1 - mf^t/j) - (k - 2k0 + 40)(1 - toO^ - (k - 4)0^(1 - H(l " "^0) 9^F 
- 4«;2 ^2(1 - to)2(^ + - e^)2a2^F - 8^2 m(l - to)01 - 0(1 - (1 - "i)0)(C + V' - ^tlj)d„,d^F 
+ Sk^ to(1 - to)(1 - toOV^I - 0^)(C + V' - ^i^)dmd^F - Ak^ 0(1 - 0^(1 ^ (1 - rn)ipfdlF 
+ 8k^ 01 - 0(1 - wO0'(1 - 0)(1 - (1 - m)i))d^d^F - Ak^ (1 - toO^0^(1 - i'fdlF . 
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Three additional equations are derived by cyclicly permuting the indices on the u variables in (6) and following 
the steps above. These three additional equations can also be obtained from (13) by the conformal symmetries of 
rectangles with arbitrary R and fixed height 1: reflection about x = R/2, reflection about y = 1/2, and reflection 
over X = y with a concurrent scaling by 1/ R in order to preserve the height. The last of these may be implemented 
by a change of aspect ratio R — > 1/R (i.e. m ^ 1 — m), exchanging x and y, and then scaling by a factor of l/R, 
so that, for example, ^ = sii{x K'\m) — > sn (y K{K' / K)\l — to) = V- The third symmetry operation introduces a 
conformal covariance factor due to the scaling, but this is absorbed into the prefactor and does not effect F . The 
three symmetry operations translate, respectively, into 

(e, V', rn) r<^ , [i. x -(\ m)!^ '"^ ' ^' ^ " "^^ " (^^^ 

Now comes a central mathematical result. There is a linear combination of the four differential equations giving 

di,diF{i,i,,Tn) = ^ . (15) 

(Since to is a variable here we have altered the usual Jacobian notation and write e.g. -F'(^, "07 ™) in place of -F'(^, ^/^Ito).) 
Thus all solutions must be of the form 

F(e,V',TO) = G(e,TO) + G(V',to) . (16) 

(An F depending on m alone is not possible, as follows from (13)). 

Now, for a given aspect ratio, ^ depends only on x and "0 depends only on y. As is spelled out in detail in [4], 
(15) then implies that the ratio of correlation functions investigated in [3] for percolation is independent of y (the 
horizontal coordinate) in the rectangle. This peculiar symmetry was, as noted in [3], first observed via simulations. 
Using (15) shows that this symmetry is exact in the continuum limit, and holds for many critical systems. 

Because of the symmetry -0, to 1 — to}, which follows from (11), it is sufficient to find the solutions of 

i^(^, ^, to) = G(^, to). The full solution set can then be obtained by letting G('(/', rn) = G('0, 1 — to) in (16). 

Inserting Gii^^m) into the differential equations and taking linear combinations allows us to cancel out all i\) 
dependence and arrive at the three equations 

= ((8 - k) (8(1 - + (8 - k)(1 - mif) + 8(6 - k)^) G - 4<to(1 - to)(k(1 - - 4(2 - 0)9mG 

+ 4/t^(l - C)(4^ -{k- 4)(1 - 2e + mi))d^G - AK^£,^m^{l - m)^dl,G (17) 

- 8K^f{l - Om(l - m)d,ndcG - AK^fil - O^d^G , 

= ((8 - k) (4 (1 - TO^^) + (8 - k)(1 - to)(1 - mO) - 8(6 - k)^(1 - 0) G 

+ 2kto(1 - to) (8 - k + 8C(1 - 2to) - 8C^ + Sm^"^ + mn^'^) d,nG (18) 
+ 2<(1 - 0(4(1 - e) + (1 - to)(-2k - 4e + 4< + (8 - K)m^))d^G 

+ 4^2^(1 - - mfd^G + 4k^^{1 - 0^(1 - to)(1 - 2$ + m£,)^„^^^G - 4k^^^{1 ~ 0^(1 - to)c'|G , and 

= ((8 - k) (8C(1 - to)(1 - ?7i^) + (8 - k)(1 - to)2) + 8(6 - k)(1 - ^f) G 

+ 4k(1 - Oto(1 - to)(4 + 4^(1 - 2to) - (1 - m)K)d„,G - 4<(1 - ^(1 - mf (4^ + k(1 - 2^))d^G (19) 

- 4^2(1 _ ^)^m^{l - mfdl^G + 8k^£,{1 - C)^m(l - mfd„,d^G - - if{l - mfdlG . 

In [3] we calculated (1), but only for the case of percolation (k = 6), and in the limit where z goes to the bottom 
edge of the rectangle. Here, because the ip dependence is entirely contained in the prefactor of (5) we expect that the 
solution space of (17)-(19) for k = 6 should contain the functions calculated in [3], up to differences in the prefactor. 
Guided by this, and with a little algebra, we find that making the substitution 

Cl _ "I'll, 3-2/11, 2 (A-m')'^/'^ 

G{i, to) = ^ ^ TT—H{m, TO = ^ ,,,, , g(TO, TO , (20) 

TO2''i,= [^(l-e)(l-mO]^ ^'-^ TO(6-«)/-[e(l-0(l-"iO]^ ' 

and taking an appropriate linear combination of the resulting equations gives the standard form of Appell's hyperge- 
ometric differential equations 

= 4) + 2[.-A)-,in-,)t^^^ ^ ^-^sd.H + .(1 - t)d.d.H , 

^ - g(s, t) - UdtH + ^^"^ ^ ~ ^ a.g + t[l - s)dtdsH + s(l - s)dlH , and (21) 

= -UtH - ^fc^a.iJ +{t- s)dtdsH , 

K K 
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with s = m and t = m^. Equations (21) have a three dimensional solution space. Among the solutions are five 
convergent Frobenius series [5]: 



/ 4 4 
Huis,t)^Fi l--;-,2 

\ K K 

( 4 4 
Hui{s,t)^Fi l--;-,2 



4 4 4 12 





>_8 


K ' 


K 


16.. 




K ' 


K 



1- sA-t 



s, t 



Hy{s,t) 



1-s, 



(s-t)i2/'--i(l_t)4/'--i / 4 4 4 12 
—Fi 1 — ; -, -; 



1 - s 



where 



s'*/k(1 - s)8/«- 



i^i (a;6i,62;c|zi,Z2) := ^ 



Ki Ki Ki t\i 



l^t' s{l-t)J ' 
(a)j+j(6i)i(62)j2i'z2-' 



i!j!(c) 



with the Pochhammer symbol (z)„ = r(z + n)/r(z), is he first of the Appell hypergeometric functions. 
This set of Appell functions allows us to write five solutions for G(^,m), valid for all k > 0, as 



Gi(C,m) 

GlIl(^,7Tl) 

Giv(C,to) 



r{2-8/K)r{m/K-l) [m(l -m)]^/''^8A-i/2 / 4 4 4 12 

-Fi 1 ; -, -; — 



r(i2/K)r(i - [(1 _ _ ^^)]4/«-i/2 

(1 - m)2/« 



Ki Hi hh Ki 



(1-m)^/'^ 

(1 - 77lOl2/'-^-3/2 

[m(l-m)]^/"-Me(l-0]^^"~'^' 
r(2-8/K)r(16/K- 1) 

r(i2/K)r(i -4/k) 



hii Hi Hi Hi 

Hi Hi Hi Hi 



1 — m, 1 — 



Hi Hi Hi Hi 



1 — TO 



1 — TO, 



1 — TO^ 



and 



(l-TO)6/«-1^4/«-l/2(l_m^)l/2 1 



m(l - f) l-( 



1 — TO^ 1 — TO^ 



(22) 
(23) 
(24) 
(25) 

(26) 



For our ranges of ^ and to values, (22)-(26) exhaust the convergent Frobenius series solutions to the differential 
equations that can be expressed with a single Fi. We can also find five other convergent Frobenius series solutions 
that can be expressed with a single Appell function of the second type, 



F2 (a;6i,&2;ci,C2|2;i,Z2) := ^ 



ia)i+jibi)i{b2)j zi"- Z2^ 
i!j!(ci)^(c2)j 



(27) 



With the definition 



n = —2 cos {Att/k) , 



(28) 
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(note that n is the parameter of the 0(n) loop models) we have 
Gvi(C: "^) ■= Gu — nGi = Giv — tiGy 

_ r(2 - 8/K)r(i6/K - i)r(4/K) 
~ r(i -4/K)r(8/K)2 

X [m(l - m)l2/« i^^ ^ 

[ 1 — 

Gvii(C,m) := Gill - nGu = Gy + (1 - n^)Gi 

r(4/K)r(2 - 8/k) [m(l - m)f/'^^s/k-i/2 



(29) 



r(8/K)r(2 - 12/k) [(1 - e)(i - mO]''/'"-!/^^ ^ 



1 



4 4 4 



hh Ki Ki Ki 



e(i-m) 

1 — 
me, 1 - e 



12 



Gviii(e, m) 



nGiii - Giv = ri(2 - n^)Gi + (r? 
r(4/K)r(l - 8/k) 

~ r(8/At- i)r(2- 12/k) 

(1 - m)2/fcmi-6A 
Gix(e, "^) := nGiii - G„ = n(2 - r?)GM + (n^ - l)Giv 



l)Gii 



.^(l-^;^,2-l^;^,2-i^ 



(30) 



(31) 



1 — TO, rn^ 



Gx(e,m) 



Gviii 
Giii - 
Gvii 



and 



1 — TO^ 

nGiv = Gi + (1 - n2)G 
1-e 



V 



(32) 



(33) 



1 — TO e 



, TO 



We derive relations (29)-(33) in appendix A. 

These Frobenius series solutions are useful, since each one is a single conformal block. It is possible to identify the 
particular block in each case by examining the leading terms, but this can be done more elegantly using the Coulomb 
gas formalism (see [4]). 

Since the solution space is three-dimensional, two of the Gs in (22)-(26) are not independent. Using the Coulomb 
gas formalism one can show that 



2Gi„ = (2 - n^)Gi + nGii + nGiv + (2 - n^)G^ 



and 



nGi + Giv = Gn + nG 



V ■ 



(34) 
(35) 



Despite this, it is convenient to consider all five solutions, as well as the alternate forms in (29)-(33) because they 
have simple interpretations in terms of physical models. 

The physical content for 0(n) loop and Q-state Potts models is examined in [4]. The normalizations of these 
solutions are chosen in part for consistency with the vertex operator formulation used there. We also show that the 
particular functions Gi, Gv, and Gvi form a natural basis of the three dimensional solution space for the critical 
models mentioned. 

The hypergeometric functions in the conformal blocks G simplify for certain k values, corresponding to various 
physical models. This is explored in [4]. 

Finally, recall that a second set of conformal block solutions that depend on -0 and to follows from Gi{ip,m) = 
G,(e,l-TO). 



B. Common functional factor and corner operators 

We next complete the transformation of the correlation function C{w) (5) from the upper half plane into the rect- 
angle using (8), by computing the common functional prefactor in the rectangle, as a function of e, ip, m (alternatively 
X, y, to) and the parameter k. 

We've chosen to write the upper half plane prefactor as 



fn{w) 



\w - u;|''i^3-2'»i/2,o [u3iUi2 



i/ii. 3-2/11.; 



{w - Ui){'W - U2)iw - U3){w - M4) 



(36) 
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This particular form was motivated by the derivation of the analogous correlation function for percolation in the 
semi-infinite strip, S = {TZ \ R — > 00} [3]. We found that it simplified the analysis to define the prefactor so that in S 

it takes the form fsiz) = [($1/2,0(^7 where ($i/2^o('2^i ^))s is the strip one point function. For percolation 

(k = 6), the exponent —11/5 = (/ii,3 — 2/i]^/2,o)/(~2/ix/2,o) is th^ ratio of leading exponents in the bulk-boundary 
fusion rules when the bulk operator approaches a free versus fixed boundary i.e. ^i/2,o{w,'w) — > w'*'^'^~^'*i/2,o0j^ 3(H) 
on a free boundary and $1/2,0(^^7'"^) ~^ v~'^'^^^^-°l{u) on a fixed boundary (see [3] or [4] for an explanation of the 
boundary conditions). Similarly, we have chosen (36) so that when mapped into the rectangle the prefactor satisfies 



fn{z) oc [($1/2,0(2, z))^ 



(/ll,3-2/li/2,o)/(-2/li/2.o) 



(37) 



The successful analysis of F(z) in the last subsection demonstrates the merit of this choice. In the remainder of this 
subsection we explicitly determine f-ji- 

Because the mapping between the upper half plane and rectangle is singular at the corners the definition of the 
boundary condition changing operators placed there needs to be adjusted. Following [6] we use the convention 



(/)^(z,) - lim(2e)- 



Uc + e 



(38) 



where <jf{zc) defines the rectangular corner operator as a boundary operator (j) approaches a corner Zc- It follows 
that the conformal weight h'^ of the rectangular corner operator is related to the weight h of the associated boundary 
operator by h'^ ~ 2h. 

Thus we may write the rectangular geometry correlation functions as 



C{Z) = {^l2iO)ri,2{R)^l2{i + RWl,2{^^l/2A^^,z))n 



lim (16£i£2e3e4) 



i(ei)'/'1.2(^ - £2)(/'l,2(i + R- £3)^1,2(1 + £4)</'l/2,o(2: z))n 



(39) 
(40) 



The prefactor in (3) is given by transforming (36) into the rectangular geometry using the transformation law 
(z) w' {z)^'(f>{w{z)) and associating the covariance factors with /t^. The result is 



fn{z) = lim 



\w'{z)\ 



2/1, 



ntik'(^K-))i 



'11,2 



-/h {w{z)) 



(16ei£2£3£4)''i-^ 

To leading order in the variables Sj we have 

{ui,U2,U3,u4 ^ {m{K'ei)^, m - m{l - m){K'e2f, 1 + {I - m){K'e3f , (K'ei)-^} 

1/2 

Using w'{z) = 2K' [w{z) (m — w{z)) (1 — w{z))] and defining 

c{m) = 2''i.^(i^')*''' '+^'''^''"(™(l - m))2''i.2 , 

we find 



(41) 



(42) 




Im w 



^Jw(ra — w){l — w) 



Im 



sn {zK'\my 



'll, 3— 2/11/2,0 



sn {zK'\m) cn {zK'\m) dn {zK'\m)\ 



We may also write f-jiiz) as 



/^(z) = c(m) 4ds(2a;i5r'|m) +4ds(2^/A"|l- 



c(m) 4(m(l-m))i/2 



i/il,3-|' '11/2,0 



t?3^ {xT:i,e~'''^) 

-TrR\ 



— |/ll.3+/ll/2,0 



Recall that /h was chosen, so that it f-ji would satisfy (37). In the upper half plane 

(4'i/2.o(«^,w))h = Imw~2''i/2,o ^ 



(43) 

(44) 
(45) 

(46) 
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and our transformation into the rectangle gives 

($1/2,0(2, - h'(z)P''l/^'°($i/2,o(w^(2),"'(2)))l 



Imif 



2K'^w {m,-w){l 



-2h 



1/2,0 



(47) 



comparing this to (43) confirms (37). In physical systems (47) can represents the density of clusters attached to a 
homogeneously wired boundary. To our knowledge no simple interpretation exists for fiz{z). 

Finally, from (44), making use of (11) and standard elliptic function properties, we find (in a slight abuse of notation) 



f{ti',m) = c{m) 



(1 - me ? 



(1 - (1 -m)V^)2 



e(i- 0(1-^0 V^(1-V)(1-(1-to)^) 



-/ll.3/2 + /li/2,0 



(48) 



This is the expression that we make use of in [4], since it uses the mathematically natural coordinates ^ and ip. 

We note that the factor in brackets in (48) is invariant under all three symmetry operations (14). Furthermore, the 
first two of these symmetries preserve m and fn{z) is completely invariant. 

The third symmetry takes to 1 — m. Combining (9), (42), and (48) we see that this implies 



(49) 



The total scaling dimension of the operators in C{z) is Ah^ ^ + 2/ix/2,05 ^i^^d thus this additional factor is the conformal 
covariance factor of a uniform scaling by So fizi^) is covariant, as required, under all symmetry operations (14). 

The factor 2^^-^ in c(to) is due to the corner operators. In general, the corner operator convention introduces 
boundary-corner fusion rules 



■3 c hS-Zio-Zii j,c 



(50) 



which modify the regular boundary-boundary fusion rules when the expansion point is a corner Xc- These corner 
OPE coefficients are related to the boundary OPE coefficients by Cf.f = 2^^02-1, where 1, 2 and 3 label the operators 
in the boundary OPE. Thus the factor of 2^^^-^ in (42) is included for consistency between the upper half-plane and 
rectangle expressions. 



III. SUMMARY AND DISCUSSION 



A. Summary 



We first summarize the main results. They include all solutions for the PDEs determining the six-point correlation 
function (1) 

C{Z) = (</)l2(O)0t2(i)*l/2,o(2,^)'^t2(^)'^t2(^ + i))K , 

in the rectangle 7^ := {z = x + i?/ G C | < a; < i?, < ?/ < 1}. Making use of the coordinates ^{x, to) and ^/'(y, m) 
from (11) 

^ = sn (x/sT'Ito)^ , and ip = sii{y K'\l — m)'^ , 
with the elliptic parameter m specifying the aspect ratio R via (10) 

(0,e-'^-«) 

777 ^ i — 

we find from (3) that any solution which involves a single conformal block can be written in the form 

C(z) = /(C,V,m)G(e,m) , 
where the algebraic prefactor / is given in (48) 



/(0,m) = c(m) 



(1 - 



(1 - (1 -to)V'')^ 



e(l-e)(l-<) 7^(l-V)(l-(l-m)V') 
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with c(m) from (42) 

c(m) = 2'''-'{K')^'''-'-+^''^/^.'>{m{l - m)f''^'^ , 

and G one of the ten solutions (22)-(26) or (29)-(33), or one of these solutions with ^ t/', m — !> 1 — m}. The 
prefactor / is independent of the physical system, while G is not. 

A bit more technically, having chosen the algebraic prefactor /e as in (5), we employ standard methods of CFT to 
give the PDEs for the conformal blocks. Because of the </>i,2S in G{z), these are second-order equations. Introducing 
the coordinates ^ and "0 (11) then leads to an interesting independence from one coordinate (15), which may be 
written 

9v.9?J^(^,0,m) = O, 

and indicates some unknown symmetry. As a consequence, conformal blocks of F only depend on two coordinates, not 
three. Some further manipulation transforms the PDEs into (21), which is the standard form for the Appell function 
Fi. A solution set, as mentioned, is given in (22)-(26) or (29)-(33). The last five of these expressions are most simply 
expressed with the Appell function F2. These solutions span a three-dimensional space; when the ip sector is included 
it is six-dimensional. The relations between the solutions are given in (34) and (35). Alternate forms for the algebraic 
prefactor / in the rectangle are given in (43)-(45). 

The results are of physical interest, since, for a variety of two-dimensional critical points, they specify the density of 
critical clusters anchored to one or both opposite ends of a rectangle with wired boundary conditions on those ends. 
These applications are discussed in [4]. 

Appendix A derives an Appell function relation that is useful in determining the Gs, and B discusses general 
conditions that must be satisfied for the conformal blocks, as here, to have a common y-dependence. 

B. Related work 

In the case k — 6, recent work by Beliaev and Izyurov [7] gives a rigorous derivation of the PDEs (7) in the case 
R ^ 00 (m — !■ 1) using SLE techniques. This provides a basis for a completely rigorous derivation of our results for 
the case of critical percolation on the triangular lattice in this limit. 
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Appendix A: Appell Function relations 

In this section we derive a relation between Appell functions (A2) used in deriving (29)-(33). The manipulations 
involved are relatively simple, involving gamma functions, Pochhammer symbols and Gauss's hypergeometric function. 
We use standard identities (see [8] for example) and assume that s,t G [0, 1). 

Beginning with the series expression for Fi we rearrange the summation with k = i + j 



where we used 

{bi)k-, _ {bi)ki-k), 
{k-j)l k\{l-bi~k)j 



Fi (a; 61, &2; c|s, si) = ^ ^ 

k=0 j=0 



{a)k{bi)k-jib2hs>^P 
{k~jy.jl{c)k 



E 

k=0 



{a)k{bi)ks'' 
kl{c)k 



3=0 



(I-61 
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which follows on rewriting the Pochhammer symbols, (z)„ = r(z + n)/T(z) and using the identity 

r(z + ?i)r(l -z-n) = (-l)'T(z)r(l -z) = (-l)"7r csc(z7r) , 

for 71 G Z. We can change the argument from i to 1 — t using standard hypergeometric function identities with the 
result that 



Fi {a;bi,b2;c\s,st) = ^ 



ia)k{bi)ks'' (1 - 6i - - k)k 
k\{c)k {l-h-k)k 



2Fi{-k,b2,bi+b2\l~t) 



k=0 
oo k 



(a)fe(&i + 62)fe(-fc),(62)js"(l-tp 



k=0 j=0 

oo oo 

EE 

oo oo 

EE 

i=0 j=o 



fc!j!(c)fc(&i + 62), 

{l+j)\jl{c),+j{bi+b2)j 

{a),+j{bi + b2),+j{b2)js'+^ {t - ly 



where the second and final expressions take advantage of the relation (1 — a — i)^ = (— l)'(a)i. 

Now we sum over the i index first. Here we use the identity (a)i+j = (a + j)i{a)j, and define the quantity 
A := c — foi — 62 — a to simplify our expressions 



Fi {a;bi,b2]c\s,st) = ^ 
3=0 



(a), (62), [sit - l]p g (a + + 62 + 



i\{c + jy 



^ (a),(b2)As(t-l)y ^ , , , 



{a),{b2)As{t-i)y ( r(-A + j)r(c + j) 



r(a + j)r(5i + 62 + i) 



(1 - s)^-^'2^^i(c - a, A + a, 1 + A - j|l - s) 



r(A-j)r(c + j) 



r(c-a)r(A + a) 

r(c)r(~A)(i^.)^ ^ (b2),(-A), 

T{a)T[bi+b2) ^^jKbi+b2)j 



Fi (1 - A - a, 1 + a - c, 1 - A + j|l - s) 



s(l™t) 



s - 1 



2i^i(c - a, A + a, 1 + A - j|l - s) 



r(c)r(A),s-- ^ {a)^{b2),{l~ty , . ^ , 

TiTi a1 2-Fi(a + J, 1 + a - c, 1 - A + j|l - 1/s) 



r(c-a)r(A + a)^ j!(l-A), 



in the final two steps we use standard identities to change the argument from s to 1 — s, and then in one case from 
1 - s to 1 - 1/s. 

If we write the functions in series form and manipulate the Pochhammer symbols we find 



r 6 ^._ r(c)r(-A)(l-.)A - - (c-a).(A + aM62),(-A) 

''^ - Tia)m+b2) E E ..,.(5,+62),(l + A-,). - ^) 



s(l-i) 



s - 1 



r(c)r(A)s-'^ 

r(c-a)r(A + a) 



EE 

j=0 1=0 



(g + j),(l + g - c),(a),(62), (1 - l/.s)^(l - ty 



r(c)r(-A)(l-s)^^(c-gMA + a), ,^ 
^ z!(l + A), ^ 



(62),(-A-z), 



r(g)r(6i+62) 

r(c)r(A)s-" 

r(c-g)r(A + g) 



s(l-0 



s - 1 



Fi{a; 1 + a - 0,62; 1 - A|l - 1/s, 1 - t) 



(Al) 
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We'll need to consider s,t € (0, 1), but neither of the terms in (Al) converge for all values in this domain. For the 
first term we perform the j summation and use the identity 



^^1 b2,-A-i,bi 



s{l-t) 



s - 1 



1 - s 
1-st 



Fi [b2,c - a + i,bi 



sil-t) 



1-st 



replacing the problematic argument with one that gives convergence in the desired range. If we then explicitly write 
the new hypergeometric series we find the first term in (Al) is 



r(c)r(-A)(i 



r(a)r(fei + 62)( 



i=0 j=0 



1 - st 



r(c)r(-A)(i - s)^+''^- f , , . , , I s{i-t) 



T{a)T{bi+b2){l- stf^- 



1 - st 



The F2 series is convergent for arguments x and y such that |a;| + |?/| < 1; it is straightforward to check that this 
expression is convergent for all s,t G (0, 1). 

Using known Fi relations we can rewrite the Appell function from the second term of (Al) as 

Fi{a;l + a-c,b2;l- A|l - l/s,l -t) = s^+''-H^+'''-''Fi{l - A - a;l + a - c,bi;l - A\l - st,l ~ t) 



which is convergent for our entire domain of {s,t}. 
Putting these pieces together we find the identity 

rfc)rf-A)fl - s)^+''2 / 
Fi (a; 61, &2; c|s, st) = , , ..^fc, c - a; A + a, 62; 1 + A, 61 + 62 



T{a)T{bi + b2){l - st) 

r{c)r{A)t'''{st)^-'- 
^ T{c-a)r{A + a) 



1 - s. 



,s(l - t) 



1-st 

Fi(l - A - a; 1 + a - c,&i; 1 - A|l - st, 1 - t) . (A2) 



Using (A2) to rewrite (22), we find Gi(C,m) = {Gu{^,m) - Gvi{^,m))/n, with Gvi(?,m) as defined in (29). The 
relations for Gvii(^,™) through Gx{(,,Tn) given in (29-33) follow similarly. 



Appendix B: Conditions for common y-dependence 



Finally we examine two related conditions for a set of correlation functions to exhibit common y dependence in the 
rectangle TZ. 

One might wonder whether, when $1/2,0 in (1) is replaced by another bulk operator with a weight 2h, the conformal 
blocks could still exhibit the ?/-independence (15). To answer this question we repeat the analysis in II A with all 
occurrences of /ii/2,o changed to h. We search for a solution F(^, tp, m) = fo{£,, ip, 'm)Fo{£,, ip, m) such that the relation 
= B^B^Fq is implied by the differential equations. The resulting conditions imply that /o must satisfy three 
differential equations, which only have a solution when h = ft.1/2.0 smd /o itself is constant. In that sense our result 
is unique. Note that the analysis depends on the second order differential equations, therefore the result is specific to 
the operators (pl 2 ■ 

One can find a necessary condition in a somewhat different way that doesn't depend on the differential equations 
and could be adapted to other corner operators. This method considers the weights of operators in the prospective 
correlation functions, along with their fusion rules. We use the operators specific to the correlation function (4), but 
this condition is easily applied to any combination of operators that one might expect to exhibit this behavior and 
for which a similar understanding of the fusion rules exists. 

Because the ratio of the correlation functions by hypothesis must not depend on the vertical position, the ratios of 
fusion products must be the same as a; — whether we approach directly from the bulk (0 < y < 1) or first move to 
a horizontal side (y = or 1). 

The argument is most easily made by referring to the 0(n) loop gas. Consider the limit a: — > 0, when the bulk point 
is taken to the left hand side of the rectangle. Then there are two cases; for one the leading order fusion product is 
the identity, for the other a four leg operator, 01,5. (This may be understood more explicitly via the treatment in 
[4].) There are, correspondingly, two possible OPEs 

0i/2,o(^, z) = Cl''"\-^''^'^-° 1 , and (Bl) 
0i/2,o(^, z) = c;(J°x'^i.-2h,/,,o . (B2) 
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We assume that this result, and the similar one below, are generally valid and not restricted to 0(n) models. 

If, on the other hand, we first bring the density operator to the bottom edge it fuses to the two- leg operator 0i_3(a;). 
Only one operator can appear in this limit for the the conformal blocks. If different operators appeared then the 
difference in scaling weights would eliminate the possibility of common y-dependence. We then move this operator 
into the lower corner of the rectangle, where there is a corner one-leg operator, 2- As mentioned in subsection II B, 
an operator in a corner of angle 7r/2 has a modified conformal weight, and thus /i^ ~ 2h. Again, from the 0(n) model 
there are two cases, either the fusion product will be to a corner one-leg operator (pl 2 or the corner three-leg operator 
(pl 4. The two possible OPEs are 

^i,3{^Wi,2iO) = Clf''^x-''^'-^^y{0) and (B3) 
0i,3(x)0t2(O) = Cl-^,^'''x''^^'--'^^-^-''^^^^(blM . (B4) 

Thus, defining p as the ratio of densities with z not allowed to belong to the left end cluster to those with z allowed 
to belong the the left end cluster, gives (depending on whether we approach via the bulk or the boundary) 

Pbuik a;''!'^ 

n, , ^2(/li,4-/ll,2) 

pbouiiclary 

Since the ratio of two correlation functions doesn't depend on which approach we take, we must have 

hi,5 = 2(/ii,4 - hi^2) (B5) 

as a necessary condition for the factorization. 

We can check this result with explicit expressions for the weights 

/ii,2 = 3/fc - 1/2 , = 15/k - 3/2 , and /ii,5 = 24/k - 2 , (B6) 

and we see that this condition is indeed true for arbitrary values of k. 



[1] A. A. Belavin, A. M. Polyakov, and A. B. Zamolodchikov. Infinite conformal symmetry in two-dimensional quantum field 

theory. Nuclear Physics B, 241:333-380, 1984. 
[2] P. DiFrancesco, P. Mathieu, and D. Senechal. Conformal Field Theory. Springer, 1999. 

[3] Jacob J. H. Simmons, Robert M. Ziff, and Peter Kleban. Factorization of percolation density correlation functions for 

clusters touching the sides of a rectangle. J. Stat. Mech: Th. Exp., 2009:P02067, 2009. 
[4] Jacob J. H. Simmons, Peter Kleban, Steven M. Flores, and Robert M. Ziff. Cluster densities at 2-d critical points in 

rectangular geometries, http://www.arxiv.com/abs/1103.5691, 2011. 
[5] P. Appell and M. J. Kampe de Feriet. Fonctions Hypergeometriques et Hyperspheriques. Gauthier-Villars, Paris, 1926. 
[6] J. L. Cardy. Conformal invariance and surface critical behavior. Nucl. Phys., B240:514-532, 1984. 
[7] D. Beliaev and K. Izyurov. A proof of a factorization formula for critical percolation. arXiv: 1011 .5822v2, 2010. 
[8] M. Abramowitz and I. A. Stegun. Handbook of Mathematical Functions. Dover Publications Inc., New York, 1964. 



